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Abstract. This is the first of a series of papers devoted to derive the metallicity of old open clusters in order 
to study the time evolution of the chemical abundance gradient in the Galactic disk. We present detailed iron 
abundances from high resolution {R ^ 40, 000) spectra of several red clump and bright giant stars in the open 
clusters IC 4651, NGC 2506 and NGC 6134. We observed 4 stars of NGC 2506, 3 stars of NGC 6134, and 5 stars 
of IC 4651 with the FEROS spectrograph at the ESO 1.5 m telescope; moreover, 3 other stars of NGC 6134 were 
observed with the UVES spectrograph on Kueyen (VLT UT2). After excluding the cool giants near the red giant 
branch tip (one in IC 4651 and one in NGC 2506), we found overall [Fe/H] values of -0.20 ± 0.01, rms = 0.02 
dex (2 stars) for NGC 2506, +0.15 ± 0.03, rms = 0.07 dex (6 stars) for NGC 6134, and -1-0.11 ± 0.01, rms = 
0.01 dex (4 stars) for IC 4651. The metal abundances derived from line analysis for each star were extensively 
checked using spectrum synthesis of about 30 to 40 Fe I lines and 6 Fe II lines. Our spectroscopic temperatures 
provide reddening values in good agreement with literature data for these clusters, strengthening the reliability 
of the adopted temperature and metallicity scale. Also, gravities from the Fe equilibrium of ionization agree quite 
well with expectations based on cluster distance moduli and evolutionary masses. 
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1. Introduction 

Our understanding of the chemical evolution of the Galaxy 
has tremendously improved in the last decade, thanks to 
the efforts and the achievements both on the observational 
and on the theoretical sides. Good chemical evolution 
models nowadays can satisfactorily reproduce the major 
observed features in the Milky Way. There are, however, 
several open questions which still need to be answered. 

One of the important longstanding questions concerns 
the evolution of the chemical abundance gradients in the 
Galactic disk. The distribution of heavy elements with 
Galactocentric distance, as derived from young objects 
like HII regions or B stars, shows a steep negative gra- 
dient in the disk of the Galaxy and of other well studied 
spirals. 

Does this slope change with time or not ? And, in case, 
does it flatten or steepen ? 

Galactic chemical evolution models do not provide a 
consistent answer: even those that are able to reproduce 
equally well the largest set of observational constraints 
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predict different scenarios for early epochs. By compar- 
ing with each other the most representative models of the 
time, Tosi (1996) showed that the predictions on the gra- 
dient evolution range from a slope initially positive which 
then becomes negative and steepens with time, to a slope 
initially negative which remains roughly constant, to a 
slope initially negative and steep which gradually flattens 
with time (see also Tosi 2000 for updated discussion and 
references) . 

From the observational point of view, the situation is 
not much clearer. Data on field stars are inconclusive, 
due to the large uncertainties affecting the older, metal 
poorer ones. Planetary nebulae (PNe) are better indica- 
tor, thanks to their brightness. PNe of type II, whose pro- 
genitors are on average 2-3 Gyr old, provide information 
on the Galactic gradient at that epoch and show gradi- 
ents similar to those derived from HII regions. However, 
the precise slope of the radial abundance distribution, and 
therefore its possible variation with time is still subject of 
debate. In fact, the PNe data of Pasquali & Perinotto 

(1993) , Maciel & Chiappini (1994) and Maciel & Koppen 

(1994) showed gradient slopes slightly flatter than those 
derived from HII regions, but Maciel, Costa & Uchida 
(2003) have recently inferred, from a larger and updated 
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PNe data set, a flattening of the oxygen gradient slope 
during the last 5-9 Gyr. 

Open clusters (OC's) probably represent the best tool 
to understand whether and how the gradient slope changes 
with time, since they have formed at all epochs and their 
ages, distances and metallicities are more safely derivable 
than for field stars. However, also the data on open clusters 
available so far arc inconclusive, as shown by Bragaglia ct 
al. (2000) using the compilation of ages, distances and 
metallicities listed by Priel (1995). By dividing her clus- 
ters in four age bins, wc find in fact no significant vari- 
ation in the gradient slope, but we do not know if this 
reflects the actual situation or the inhomogeneities in the 
data treatment of clusters taken from different literature 
sources. Inhomogeneity may lead, indeed, to large uncer- 
tainties not only on the derived values of the examined 
quantities, but also on their ranking. Past efforts to im- 
prove the homogeneity in the derivation of abundances, 
distances or ages of a large sample of OC's include e.g., 
the valuable works by Twarog et al. (1997) and Carraro 
et al. (1998), but in both cases they had to rely on liter- 
ature data of uneven quality. The next necessary step is 
to collect data acquired and analyzed in a homogeneous 
way. 




Fig. 1. Field of 10 x 10 arcmin^ centered on NGC 2506, 
with the 4 stars observed with FEROS indicated by their 
numbers according to Marconi et al. (1997) 



For this reason, we are undertaking a long term project 
of accurate photometry and high resolution spectroscopy 
to derive homogeneously ages, distances, reddening and 
element abundances in open clusters of various ages and 
galactic locations and eventually infer from them the gra- 
dient evolution. Age, distance, and reddening arc ob- 
tained deriving the colour-magnitude diagrams (CMDs) 
from deep, accurate photometry and applying to them the 
synthetic CMD method by Tosi et al. (1991). Accurate 
chemical abundances are obtained from high resolution 
spectroscopy, applying to all clusters the same method of 
analysis and the same metallicity scale. 

Up to now wc have acquired the photometry of 25 clus- 
ters and published the results for 14 of them (sec Bragaglia 
et al. 2004 and references therein, Kalirai & Tosi 2004). We 
have also obtained spectra for about 15 OC's, and have 
already obtained preliminary results for a few of them. 
Complete spectroscopic abundance analysis has been pub- 
lished so far only for NGC 6819 (BragagUa et al. 2001). 
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In this paper we present our results for three OC's, 
namely NGC 2506, NGC 6134, and IC 4651. Observational 
data are presented in Section 2, while Section 3 is devoted 
to the derivation of atmospheric parameters, equivalent 
widths and iron abundances. A check of the validity of our 
temperature scale is derived from comparison with red- 
dening estimates from photometry in Section 4. Spectral 
synthesis for all stars and its importance to confirm the 
validity of our findings is discussed in Section 5. Finally, 
Sections 6 and 7 present a comparison with literature de- 
terminations, and a summary. 



Fig. 2. Field of 10 x 10 arcmin^ centered on NGC 6134, 
with the 6 stars observed with FEROS or UVES indicated 
by their numbers according to Lindoff (1972). 

2. Observations 

Stars were selected on the basis of their evolutionary 
phases derived from the photometric data. We mainly tar- 
geted Red Clump stars (i.e. stars in the core Helium burn- 
ing phase) because, among the evolved population, they 
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Fig. 3. Field of 10 x 10 arcmin^ centered on IC 4651, 
with the 5 stars observed with FEROS indicated by their 
numbers according to Lindoff (1972), and Meibom (2000). 



are the most homogeneous group, and their temperatures 
are in general high enough that, even at the high metal- 
licity expected for open clusters, model atmospheres can 
well reproduce the real atmospheres (this seems to be not 
strictly valid for stars near the Red Giant Branch - RGB 
- tip). We considered only objects for which membership 
information was available (except for one object in NGC 
6134, which however turned out to be member a poste- 
riori), from astrometry (NGC 2506, Chiu & van Altena 
1981), and/or radial velocity (NGC 2506: Friel & Janes 
1993, Minniti 1995; NGC 6134: Claria & MermiUiod 1992; 
IC 4651: Mermilliod et al. 1995). 

The three clusters were mainly observed with the 
spectrograph FEROS (Fiber-fed Extended Range Optical 
Spectrograph) mounted at the 1.5m telescope in La Silla 
(Chile) from April 28 to May 1 2000. FEROS is bench 
mounted, and fed by two fibers (object -f- sky, entrance 
aperture of 2.7 arcsec). The resolving power is 48000 and 
the wavelength range is AA 370-860 nm. There is the pos- 
sibility of reducing the data at the telescope using the 
dedicated on-line data reduction package in the MIDAS 
environment, so we could immediately obtain wavelength 
calibrated spectra. Multiple exposures for the same star 
were summed. 

Three additional stars in NGC 6134 were observed in 
July 2000 with UVES (UV- Visual Echelle Spectrograph) 
on Unit 2 of the VLT ESO-Paranal telescope, as a backup 
during the execution of another programme (169.D-0473). 
These data were acquired using the dichroic beamsplitter 
#2, with the CD2 centered at 420 nm (spectral cover- 
age AA 356-484 nm) and the CD4 centered at 750 nm 



(AA 555-946 nm). The slit length was 8 arcsec, and the 
slit width 1 arcsec (resolution of 43000 at the order cen- 
ters). The UVES data were reduced with the standard 
pipeline which produces extracted, wavelength calibrated, 
and merged spectra. 

Finding charts for all observed targets are in Figure ^ 
Figure|21and Figure|21 The evolutionary status of observed 
stars is indicated by their position in the CMD, as shown 
in Figure 0] 

Table ^ gives a log of the observations, and more in- 
formation on the selected stars is listed in Table El where 
values of the S/N are measured at about 670 nm (for mul- 
tiple exposures, these values refer to the final, co-added 
spectra). 

3. Atmospheric parameters and iron abundances 

3.1. Atmospheric parameters 

We derived effective temperatures from our spectra, by 
minimizing the slope of the abundances from neutral Fe I 
lines with respect to the excitation potential. The gravities 
(logg) were derived from the iron ionization equilibrium; 
to this purpose, we adjusted the gravity value for each 
star in order to obtain an abundance from singly ionized 
lines of iron 0.05 dex lower than the abundance from neu- 
tral lines. This was done to take into account the same 
difference between Fe II and F I in the reference analy- 
sis made using the solar model from the Kurucz (1995) 
grid, with overshooting: the present adopted values are 
logn(Fe)=7.54 for neutral iron and 7.49 for singly ionized 
iron. 

As usual, the overall model mctallicity [A/H] was cho- 
sen as that of the model atmosphere (extracted from 
the grid of ATLAS models with the overshooting option 
switched on) that best reproduces the measured equiva- 
lent widths (EW's). 

Finally, we determined the microturbulent velocities 
assuming a relation between log g and Vt . This was found 
to give more stable values than simply adopting individual 
values of Vt by eliminating trends in the abundances of 
Fe I with expected line strengths. In fact, we found that 
the star-to-star scatter in Fe abundances was reduced by 
adopting such a relationship. 

Our adopted atmospheric parameters for the three 
clusters are listed in Tabled 

3.2. Equivalent widths 

We measured the EW^s on the spectra employing an 
updated version of the spectrum analysis package devel- 
oped in Padova and partially described in Bragaglia et 
al. (2001) and Carretta et al. (2002). Measurements of Fe 
lines were restricted to the spectral range 5500-7000 A to 
minimize problems of line crowding and difficult contin- 
uum tracing blueward of this region, and of contamination 
by telluric lines and possible fringing effects redward. 
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Table 1. Log of the observations. ID is taken from Marconi et al. (1997) for NGC 2506, from Bruntt et al. (1999) for 
NGC 6134, and from Lindoff (1972) for IC 4651. IDbda is the identification number used in the BDA (Mermilliod 

1995). Coordinates are at J2000, and exposure time is in seconds. F indicates FEROS spectra, and U indicates UVES 
spectra. The P values given for the stars in NGC 2506 arc membership probabilities. 
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Table 2. Data for the observed stars in NGC 2506, NGC 6134, and IC 4651. V, B-V are: CCD magnitudes from 
Marconi et al. (1997) for NGC 2506; photoelectric measurements from Claria & MermilUod (1992) for NGC 6134; 
photoelectric measurements from Eggcn (1971), Lindoff (1972), and Jennens & Heifer (1975), as reported in the BDA, 
for IC 4651. K is taken from 2MASS (Cutri et al. 2003). Stromgren y, b-y are: from Bruntt et al. (1999) for NGC 
6134, and Anthony- Twarog & Twarog (2000) for IC 4651. S/N has been measured at about 670 nm. Radial velocities 
(column 8, in km s~^), are heliocentric. 
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Fig. 4. CMD's for the three clusters, with the stars analyzed here indicated by larger symbols, (a) NGC 2506: B and V 

taken from Marconi et al. (1997); only stars with membership probability greater than 0.7 (Chiu & van Altcna 1981), 
are shown, (b) NGC 6134: V, b-y taken from Bruntt et al. (1999). (c) IC 4651: V, b-y taken from Anthony-Twarog & 
Twarog (2000). 

Table 3. Adopted atmospheric parameters and derived Iron abundances for observed stars in NGC 2506, NGC 6134, 
and IC 4651; n indicated the number of lines retained in the analysis. In the last column we give for comparison the 
value of gravity derived from photometric information. 
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At these high mctallicitics and cool temperatures, con- 
tinuum tracing is a major source of uncertainty at the res- 
olution of our spectra, hence we chose to use an iterative 
procedure. For all clusters, the fraction of the 200 spec- 
tral points centered on each line to be measured and used 
to derive the local continuum level was first set to 1/4. 
Following an extensive comparison with spectrum synthe- 
sis of Fe I and Fe II lines (see below) , we found that the 
abundances for the stars in NGC 2506 and NGC 6134 
(from FEROS spectra) were much too large. This was ex- 
plained by a too high continuum tracing in these spectra; 
for these two sets of data we finally adopted 1/2 as the 



fraction of spectral points to be used in the estimate of the 
local continuum. With this new parameters, we performed 
again the automatic measurements of EW's. Then, the 
procedure strictly followed that described in Bragaglia et 
al. (2001). 

We used stars in the same evolutionary phase (in 
this case, clump stars) to obtain an empirical estimate 
of internal errors in the EW's. We performed the cross- 
comparisons of the sets of EW's measured for clump stars: 
3 stars in IC 4651, 2 in NGC 2506 and 6 in NGC 6134 (in 
the latter case we treated separately the 3 stars observed 
with FEROS and the 3 from UVES). Assuming that er- 
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rors can be equally attributed to both stars in the couple 
under consideration, we obtain typical errors in EW^s of 

2.5 mA for clump stars in IC 4651, 2.7 mA in NGC 2506, 

3.6 mA and 3.1 mA in NGC 6134 (FEROS and UVES 
spectra, respectively) . 

The classical formula by Cayrel (1989) predicts that, 
given the full width half-maximum and the typical S/N 
ratios (see Table 12)1 of our spectra, the expected errors are 
1.5, 1.8, 2.1 and 0.8 mA, respectively, in the four cases. 

The comparison with the observed errors shows that 
another source of error (quadratic sum) has to be taken 
into account. This is likely to be uncertainties in the 
positioning of the continuum, an ingredient neglected in 
Cayrel's formula. If we consider for simplicity's sake lines 
of triangular shape and the relationship between FWHM 
and central depth used by our automatic procedure, we 
can estimate that the residual discrepancy can be well ex- 
plained by errors in the automatic continuum tracing at 
a level of 1% for clump stars in all clusters. 

Sources of oscillator strengths and atomic parameters 
are the same of Gratton et al. (2003) and discussion and 
references are given in that paper. 

3.2.1. Estimate of errors in atmospheric parameters 

Errors in effective temperatures: Uncertainties in the 
adopted spectroscopic temperatures can be evaluated 
from the errors in the slope of the relationship between 
abundances of Fe I and excitation potentials of the lines. 
Varying by a given amount the effective temperature, 
these errors allow us to estimate a standard error of 88 ±9, 
rms = 29 K (from 11 stars, excluding stars at the RGB 
tip). Hence, we adopt a standard error of about 90 K, 
corresponding to an average rms of 0.018 dex/eV in the 
slope. 

This error {a■^^^ — errand ~^ "'syst.) each individual 
star includes two contributions: 

— a random term cr^^^^ , reflecting the various sources of 
errors (Poisson statistics, read-out noise, dark current) 
of spectra, that affects the measurements of the EW 

— a term crgygt due to a number of effects that are sys- 
tematic for each given line (e.g. the presence of blends, 
of possible features not well accounted for in the spec- 
tral regions used to determine the continuum level, the 
line oscillator strength, etc.) 

The first appears as a star-to-star scatter both in the 
slope abundance/excitation potential and in the error as- 
sociated to this slope. Analogously, the second contribu- 
tion shows up as a systematic uncertainty both in the 
slope and in the associated error. In order to estimate the 
true random internal error in the derived temperatures, 
this systematic contribution has to be estimated and sub- 
tracted. 

We proceed in the following way. Let j be the index 
associated to the stars and N the number of stars; let i be 
the index associated to the lines and Mj the total number 
of lines of neutral iron used in the j— th star. We then 



have that the average variance of the distribution of the 
abundances from individual lines for each star is: 

j Mpr 



and as an estimate of the systematic contribution we 
may take: 

Nj-l 



In our case, if we exclude from the computations the 
two tip stars and star 456 in NGC 2506 (the one with lower 
S/N), we obtain that (Xtot = 0.141 dex and Csyst. = 0.110 
dex. Thus the truly random contribution to the standard 
deviation is (quadratic sum) 0.088 dex, hence the frac- 
tion of the standard deviation due to star-to-star errors is 
0.088/0.141= 0.624. Hence, we may expect that the inter- 
nal, random error in the temperatures is actually ~ 55 K, 
to be compared with the observed value as derived from 
independent methods (see below). 

Errors in surface gravities: Since our derivation of atmo- 
spheric parameters is fully spectroscopic, a source of in- 
ternal error in the adopted gravity comes from the total 
uncertainty in Tcff. 

To evaluate the sensitivity of the derived abundances 
to variations in the adopted atmospheric parameters for 
Fe (reported in Table 21 , we re- iterated the analysis of 
the clump star 146 of IC 4651 while varying each time 
only one of the parameters of the amount corresponding 
to the typical error, as estimated above. Considering the 
variation in the ionization equilibrium given by a change 
of 90 K, aA/9Toff = 0.143 dex, and by a change of 0.2 
dex, dA/dg — 0.096 dex, an error of 90 K in tempera- 
tures translates into a contribution of 0.298 dex of error 
in gravity. 

From the above discussion, the contribution due to the 
random error in temperature and to the error in gravities 
is then of 0.18 dex. 

A second contribution is due to the errors in the mea- 
surements of individual lines. If we assume as the error 
from a single line the average rms of the abundances from 
Fe I lines (0.14 dex, considering the same stars as above), 
again we must take into account only the random contri- 
bution, 0.088 dex, as estimated above. Weighting this con- 
tribution by the average number of measured line (n=lll 
for Fe I and 12 for Fe II) we obtain a random error of 0.03 
dex for the difference between the abundance from Fe II 
and Fe I lines. This corresponds to another 0.06 dex of 
uncertainty in log g, to be added in quadrature to the pre- 
vious 0.18 dex. The total random internal uncertainty in 
the adopted gravity (quadratic sum) is then 0.20 dex that 
compares well with the estimate derived from the method 
described in the following. 
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Errors in microturbulent velocities: To estimate the proper 
error bars in the microturbulent velocity values, we started 
from the original atmospheric parameters adopted for the 
clump star 146 in IC 4651. In this star the slope of the line 
strength vs abundance relation has the same value of the 
quadratic mean of errors in the slope of all stars, Hence it 
may be considered as a typical error. The same set of Fe 
lines was used to repeat the analysis changing vt until the 
IfT value from the slope of the abundance/line strength re- 
lation was reached. A simple comparison allows us to give 
an estimate of la errors associated to Vt'. they are about 
0.1 km s~^, thanks to the large number of used lines of 
different strengths. Since we estimated that random errors 
are about 62% of the error budget, the random internal 
error is then 0.06 km s^^. 

In order to estimate the total error, we have to recall 
that our final Vt values are derived using a relationship 
between gravities and microturbulent velocities from the 
former full spectroscopic analysis: Vt = 1.5 — 0.13 x \ogg. 
The rms scatter around this relationship is 0.17 km s^^, 
so that the final error bar in microturbulent velocities de- 
rived through the relation is 0.16 km s^^, which includes 
a (small) component due to errors in gravity and another 
component of physical scatter intrinsic to the different 
stars. 

Col. 7 of Table 01 allows to estimate the effect of errors 
in the EW; this was obtained by weighting the average 
error from a single line with the square root of the mean 
number of lines (listed in Col. 6 of the Table) measured 
for each element. 

Notice that Cols. 2 to 5 in Table El are only meant to 
evaluate the sensitivity of the abundances to changes in 
the adopted atmospheric parameters. As discussed above, 
the actual random errors involved in the present analysis 
are, more reasonably, 50 K in Toff, 0.2 dex in logg, 0.05 
dex in [A/H] and 0.16 km s^^ in the microturbulent veloc- 
ity. Moreover, the sensitivities in Cols. 2 to 5 of this Table 
are computed assuming a zero covariance between the ef- 
fects of errors in the atmospheric parameters. In principle, 
this assumption is not strictly valid, since there are cor- 
relations between different parameters: (i) effective tem- 
perature and gravities are strictly correlated, so that for 
each 90 K change in Toff there is a corresponding change 
of 0.3 dex in logg, since gravities are derived from the 
ionization equilibrium of Fc; (ii) there is a correlation be- 
tween Vt and Toff, since lines of low excitation potential 
tend to be systematically stronger than lines with high 
E.P. Thus, a change in vt gives a change in effective tem- 
peratures as derived from the excitation equilibrium. The 
sensitivity of this effect is not very large, anyway, and it 
is about 20 K for each 0.1 km s~^ change in Vf. 

Taking into account these correlations in the compu- 
tation of the total uncertainty in the abundances, we have 
an error of 0.040 dex due to the uncertainty of 0.16 km s~^ 
in the vt value; moreover, another uncertainty of 0.032 dex 
comes from the uncertainty in temperature. If we neglect 
the other contributions (that are small, anyway) from the 
quadratic sum we derive that the abundance of an indi- 



vidual star bears a random error of 0.051 dex (listed in 
the last Col. of Table O that compares very well with the 
observed scatter (0.053 dex) of individual stars around the 
mean value observed for each single cluster. 

In the Table we omit in the last Column the total 
random error in the Fe II abundance, since we force it 
to be identical to that of Fe I. 

3.2.2. Gravities from stellar models 

We can compare the gravities derived solely from the spec- 
troscopic analysis with what we obtain from the photo- 
metric information. We can compute gravities using the 
relation: log g = -10.607 -I- io5(M/M0) - Log^L/Lo) + 
4 X Log(Toff). For the temperatures, that enter the rela- 
tion both directly and indirectly through the bolometric 
correction (BC), we use our derived ones, since they are on 
the Alonso et al. (1996) scale, and are therefore very close 
to what would be obtained from dereddened colours (see 
next section). Values for masses are obtained reading the 
Turn-Off values on the Girardi et al. (2000) isochrones 
for solar metallicity at the generally accepted ages for 
these clusters (1.7 Gyr and M=1.69 M© for NGC 2506 
and IC 4651, 0.7 Gyr and M=2.34 Mq for NGC 6134). 
Adoption of any reasonable different age or isochrone, and 
the fact that we are dealing with (slightly) more massive 
stars, since they have already evolved from the main se- 
quence, would have a negligible impact on the final gravity. 
Absolute magnitudes are computed using the literature 
distance moduli and reddenings: (m-M)o = 12.6 (Marconi 
et al. 1997), 9.62 (Bruntt et al. 1999), and 10.15 (Anthony- 
Twarog & Twarog 2000) for NGC 2506, NGC 6134, and 
IC 4651 respectively. The BC is derived from eqs. 17 and 
18 in Alonso et al. (1999), and we assume Mf,oi,0 = 4.75. 
Results of these computations are presented in last column 
of Table O 

If we consider only the red clump stars, the average 
difference between the gravities derived from photometry 
and spectroscopy for the whole sample is 0.02 ±0.07 [rms 
— 0.24, 12 stars). When we do the same for each cluster, 
we find instead: +Q.21 ±0.07 (rms = 0.09, 2 stars) for 
NGC 2506, -0.16 ±0.07 {rms = 0.17, 6 stars) for NGC 
6134, and +0.18 ±0.11 (rms = 0.22, 4 stars) for IC 4651. 
The internal error for each cluster 0.09 dex) is per- 
fectly compatible with the internal errors in temperature 
and Fe abundances. Instead, the cluster to cluster disper- 
sion (0.21 dex) appears larger. Many factors - and their 
combinations - could contribute to this dispersion, among 
these: internal errors (0.09 dex), errors in distance mod- 
uli (a conservative estimate of about 0.2 mag translates 
into about 0.08 dex) and reddenings (almost negligible, 
since the values tied to the adopted distance moduli are 
very similar to ours), in ages (hence masses, giving a small 
contribution of less than about 0.03 dex even for a 20 % 
age error), non homogeneity of the data sources, small dif- 
ferences in the helium content of these clusters not taken 
into account in the distance/age derivation (responsible 
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Table 4. Sensitivities of abundance ratios to errors in the atmospheric parameters and in the equivalent widths 



Ratio 


(+90 K) 


A log 3 
(+0.2 dex) 


A [A/H] Avt 
(+0.1 dcx) (+0.1 km s-1) 


<N > 


A EW 


tot. 

(dex) 








Star IC4651-146 (Clump) 








[Fe/H]I 
[Fe/H]II 


+0.057 
-0.086 


+0.011 
+0.107 


+0.010 -0.044 
+0.035 -0.043 


111 
12 


+0.013 
+0.040 


0.051 



for less than about 0.05 dex in logg), systematic effects 
somewhat tied to age (the two clusters with similar ages 
have also similar differences in gravities, and opposite in 
sign to the younger one), etc. Moreover, plotting the de- 
rived values of photometric gravities vs our spectroscopic 
ones (Fig. [S)), we see a trend with temperature: clump 
stars of different clusters are segregated due to differences 
in metallicities and ages. There is then a hint that the sys- 
tematic effect affecting the spectroscopic gravities depends 
on the clump temperature. What is the real physical effect 
is not clear: it might be due to some Fe 11 lines blended 
at lower temperature but clean at higher temperatures, or 
to the atmospheric structure of stars systematically vary- 
ing with temperature. Recently, various authors (see e.g. 
Allende-Prieto et al. 2004, and references therein) showed 
that similar problems are also present in the analysis of 
dwarfs; however, in our program stars the effect is much 
more smaller and it does not significantly affect our con- 
clusions on reddening, temperature and metallicity scales. 

Given the small size of our sample, we cannot disen- 
tangle the above possibilities, and we postpone a more in 
depth discussion until we have examined more clusters in 
our sample. 

On the other hand, for stars at the RGB tip the 
agreement between the gravities computed with the two 
methods is much worse, with gravities derived from spec- 
troscopy lower than those obtained from evolutionary 
masses. 

In summary, when considering only clump stars, the 
agreement between spectroscopic gravities (i.e. obtained 
through the ionization equilibrium) and evolutionary 
gravities is on average good. This result leads to two rel- 
evant conclusions: (i) if our gravity values are correct, 
also our temperatures must be correct, hence our find- 
ings strongly support our temperature scale, which in turn 
translates in a reliable metallicity scale; (ii) the small dif- 
ferences in gravity found with the two methods seem to 
hint convincingly against the existence of large departures 
from the LTE assumption (see Gratton et al. 1999): high 
quality spectra coupled with reliable oscillator strengths 
from laboratory experiments are able to provide solid Fe 
abundances for moderately warm stars. 

As far as the cooler RGB tip stars are concerned, the 
situation is clearly more complex and less stable. In this 
case, we cannot exclude departures (even relevant) from 
LTE, since the lower atmospheric densities do favour non- 
local effects. Moreover, it is possible that the atmospheric 
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Fig. 5. Spectroscopic gravities as a function of the derived 
photometric gravities. Different symbols are used for the 
3 clusters. In the inset, only the clump stars and red giant 
stars are plotted, excluding the stars near the RGB-tip. 

structure of giant stars near the RGB tip is not well re- 
produced by presently existing model atmospheres of the 
Kurucz grid (see Dalle Ore 1993). For these reasons we be- 
lieve it is safer to concentrate on warmer objects for this 
kind of analysis, in order to obtain more reliable results. 
In old OC's the red clump stars are the best choice, given 
their temperatures and luminosities. 

4. Reddening estimates from spectroscopy 

Our effective temperatures are derived entirely from our 
spectra. As such, they are reddening-free and this ap- 
proach allows us to give an estimate of the reddening to- 
ward the clusters which is independent of the photometric 
determinations. For this exercise, we used only the stars 
with Tefi < 4400 K, i.e. on the clump or near it. 

First we collected all the available broad and 
intermediate-band photometry for the analyzed stars. 
Johnson and Stromgren photometry was taken from the 
references indicated in Table 2, while JHK photometry 
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Table 5. Values for the reddening as determined in the 
present paper; the three E(B-V) values for each cluster 
derived from individual colours, i.e. from (B-V), (b-y), and 
(V-K) respectively, are also given. The last line gives the 
adopted reddening. 





NGC 2506 


NGC 6134 


IC 4651 


[Fe/H] 


-0.20 


+0.15 


+0.11 


Teft(ref.) 


5000 K 


5000 K 


4500 K 


E(B-V)i3v 


0.058±0.016 


0.355±0.005 


0.087+0.022 


E(B-V)6, 




0.340±0.006 


0.082+0.015 


E(B-V)vK 


0.076±0.007 


0.388±0.005 


0.080 


E(B-V) adop. 


0.073±0.009 


0.363±0.014 


0.083+0.011 



was obtained from the 2MASS survey (from the Point 
Source Catalogue of the All-Sky Data Release, found 
at http://www.ipac.caltech.edu/2mass/) and trans- 
formed to the TCS system. 

We then adopted the colour-temperature transforma- 
tions by Alonso et al. (1996) and we entered our spectro- 
scopic Toff values in these relations to obtain de-reddened 
colours. Comparison with the observed colours provides 
an estimate of the E{B — V) value. 

Our final reddening values are derived as the weighted 
average of the reddening values as given from individ- 
ual colours (adopting E{b - y) = 0.72E{B - V) and 
E{V -K) = 2.75E{B - V): Cardelh et al. 1989). As re- 
sulting errors, we adopted the larger between the internal 
error and the spread in the values obtained from individ- 
ual colours. 

Our results for the three clusters are listed in Table [3 
The agreement with previous literature data is quite good 
(see Tableinil. Overall, the differences are about 0.01 mag, 
with our values being slightly lower, on average by 0.008 
mag. This implies that our temperatures might be a little 
too low: 0.01 mag in B — V translates into an offset of 
about 20 K in the derived temperatures. 

Our temperatures for individual stars have attached 
random errors of about 40-50 K; since we use from 3 to 
6 stars in each cluster, we expect errors of about 20-30 
K per cluster, on average, hence 0.010-0.015 mag in the 
E{B — V) value. Notice that the ~ 90 K estimated as 
the total internal error in temperature (see above, Sect. 
3.2.1) do include not only a random component due to 
measurement errors, but also a systematic component due 
to the set of lines used (e.g. some lines always giving too 
high or too low abundances). 

On the other hand, such a good agreement with liter- 
ature data strengthens the reliability of our temperature 
scale, which is only 15 K (with 10-15 K of uncertainty) 
too low with respect to that by Alonso et al. (1996) from 
IRFM. This, in turn, might translate in our abundances 
being underestimated by (only) 0.01 dex. 



Table 6. Literature values for E(B-V) and E(b-y). 
SFD98=Schlegel et al. 1998; M+97=Marconi et al. 1997; 
K+01=Kim et al. 2001; D+02=Dias et al. 2002 (cata- 
logue); KF91=Kjeldsen & Frandsen 1991; CM92=Claria 
& Mermilliod 1992; B+99=Bruntt et al. 1999; 
AT+88=Anthony-Twarog et al. 1988; ATTOO=Anthony- 
Twarog & Twarog 2000; M+02=Meibom et al. 2002 



E(B-V) 


E(b-y) 


Reference & notes 






NGC 2506 




0.087 




SFD98 




0-0.07 




M+97 [Fe/H] ~ 


to - 0.4 


0.04+0.03 




K+01 




0.081 




D+02 [Fe/H] = 
NGC 6134 


-0.37 


0.46+0.03 




KF91 




0.35+0.02 




CM92 [Fe/H] = 


-0.05 


0.365+0.006 


0.263+0.004 B+99 [Fe/H] = 


+0.28 


0.395 




D+02 [Fe/H] = 
IC 4651 


+0.18 


0.086 




AT+88 




0.241 




SFD98 




0.086 


0.062 


ATTOO [Fe/H] = 


= +0.08 


0.099 


0.071 


ATTOO [Fe/H] = 


= +0.12 




0.071 


M+02 [Fe/H] ~ 


+0.12 



5. Analysis with synthetic spectra: a check 

When dealing with high-metallicity, rather cool stars, the 
line crowding might become rather difficult to treat simply 
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Fig. 6. Comparison between the observed spectrum of the 
field metal-rich giant star HR 3627 (open squares) and the 
synthetic spectrum computed by optimizing the Kurucz 
line lists (solid lines). Small regions around 2 Fe I and 2 
Fe II lines are shown. 
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with the standard method of Une analysis. In order to 
check our results, we devised a new procedure based on 
an extensive comparison with spectrum synthesis: 

— as a first step, line lists were prepared for selected re- 
gions of 2 A centered on 44 Fe I lines and 9 Fe II lines 
chosen among those employed in our EW analysis. 
These lines are in the wavelength range between 5500 
and 7000 A, and were selected because they have very 
precise log (;/'s, and are in comparatively less crowded 
spectral regions. 

— For each region around Fe lines we computed a syn- 
thetic spectrum, starting from the Kurucz line list. 
Our list was then optimized by comparing it with 
the very high resolution (i? ^ 140,000), high signal- 
to-noise spectrum of star HR 3627, a metal rich 
([Fe/H]=-f 0.35) field giant. This star was selected be- 
cause its spectrum is very line rich, due to the com- 
bination of low temperature and high metal abun- 
dance. Hence we expect that a line list compiled from 
this spectrum would be adequate to discuss even the 
most difficult cases among the stars in our sample. 
HR 3627 was observed with the SARG spectrograph 
mounted at the TNG, and its atmospheric parameters 
were derived from a full spectroscopic analysis and are: 
4260/1. 95/-H0.35/1. 20 (effective temperature, gravity, 
metal abundance and microturbulent velocity, respec- 
tively) . An example of the optimized synthetic spectra 
for star HR 3627 is given in Figure where two Fe I 
lines and two Fe II lines are showed, with the corre- 
sponding synthetic spectrum overimposed^. 

— After the line lists were optimized to well reproduce the 
spectrum of star HR 3627 (in particular for contami- 
nants like CN, with large contribution in cool, metal- 
rich stars), we ran an automatic procedure that allows 
a line-by-line comparison of the observed spectrum and 
the synthetic spectra of the Fe lines. 

— For each program star, 7 synthetic spectra were com- 
puted for each Fe line, varying the [Fe/H] ratios from 
—0.6 to -1-0.6 dex, in step of 0.2 dex. 

— The EW's of these synthetic spectra were measured 
by using a local pseudo-continuum that considers 2 
small regions 0.8 A wide on each side of the line. The 
measurements are then saved. 

— The same lines are measured exactly in the same way 
on the observed spectra of the program stars. 

— A line-to-line comparison is made between the ob- 
served EW^s and those measured on the synthetic 
spectra; the synthetic EW that best matches the ob- 
served one is then adopted as the one giving the right 
abundance; examples of the matches obtained for star 
421 in NGC 6134 are shown in Figure 0. 

— Finally, the resulting set of abundances are treated 
analogously to the traditional line analysis, tempera- 
tures and microturbulent velocities are derived by min- 

^ Optimized line lists are available upon request from the 
first author 




Fig. 7. Comparison between the observed UVES spec- 
trum of star 421 in NGC 6134 (open squares) with syn- 
thetic spectra of two Fe I and two Fe II lines (solid lines). 
The spectral synthesis was computed by using the atmo- 
spheric parameters appropriate to the star and 3 Fe abun- 
dances: [Fe/H] =-0.2, 0.0 and -f 0.2 dex from top to bottom, 
respectively. 

imizing trends of abundances as a function of the ex- 
citation potential and line strengths and so on. 

As mentioned above, the first iteration using this tech- 
nique showed a discrepancy with the values derived by 
line analysis for the NGC 2506 stars. For stars in NGC 
6134 the disagreement was found for stars with FEROS, 
but not with UVES spectra, where the adopted fraction 
of high spectral points around the lines used for the lo- 
cal continuum tracing was higher (1/2 rather than 1/4). 
This result leads us to conclude that the culprit was the 
location of the continuum, traced too high in the FEROS 
spectra of NGC 6134 and NGC 2506. We then repeated 
the EW measurements by adopting 1/2 as the fraction in 
the routine. This time, the improvement was significant 
and the agreement is fairly good, as shown in Table [T] 

The larger scatter in the Fe I abundances as obtained 
from spectrum synthesis is likely due to the method used 
to measure the local pseudo continuum around the Fe 
lines. 

6. Results and discussion 

Using the values derived for [Fc/H] only for the clump 
stars (Table ^ we have the following iron abundances: 
[Fc/H] = -0.20 (ct= 0.02, 2 stars) for NGC 2506, -f0.15 
(cr= 0.07, 6 stars) for NGC 6134, and +0.11 {a= 0.01, 3 
stars) for IC 4651. 
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Table 8. Literature values for metallicity, and methods used. The two values given for NGC 2506 by Marconi et al. 
correspond to the best reproductions of the observed CMDs with synthetic ones based on the Padova tracks (Bressan 
et al. 1993) at Z=0.008 and 0.02; note that further unpublished analyses with updated evolutionary tracks have 
shown that the solar solution should be excluded. The two solutions for IC 4651 given by Anthony-Twarog & Twarog 
correspond to two different calibrations of the intrinsic b-y colour versus metallicity. 



Reference 


[Fe/H] 


E(B-V) 


method 




NGC 2506 






this paper 


—0.20 


0.073 


high-res sp. 


Geisler et al. 1992 


-0.58 




Washington ph. 


Friel & Janes 1993 


—0.52 


0.05 


low-res sp. 


Marconi et al. 1997 


~ —0.4 to 0.0 


0.05 to 


synthetic CMD 


Twarog et al. 1997 


—0.38 


0.05 


low-res sp.+ DDO 


Iriel et al. 2002 


—0.44 


0.05 


low-res sp. 










this paper 


+0.15 


0.363 


high-res sp. 


Claria & Mermilliod 1992 


-0.05 


0.35 


Wash. + DDO 


Twarog et al. 1997 


+0.18 


0.35 -0.39 


DDO 


Bruntt et al. 1999 


+0.28 


0.365 


Stromgren 




IC 4651 






this paper 


+0.11 


0.083 


high-res sp. 


Twarog et al. 1997 


+0.09 


0.11 - 0.12 


DDO 


Anthony-Twarog & Twarog 2000 


+0.077 


0.086 


Stromgren 




+0.115 


0.099 


Stromgren 


Meibom et al. 2002 


~ +0.12 


0.10 


CMD (Yale) 



Table 7. Comparison between Fe I abundances from spec- 
trum synthesis and EW analysis for observed stars in 
NGC 2506, NGC 6134, and IC 4651. 



Star 


[Fe/H]I 


rms 


n 


[Fe/H] I 


Phase 




SS 






EW 








NGC 


2506 






459 


-0.42 


0.264 


30 


-0.37 


RGBtip 


438 


-0.21 


0.144 


26 


-0.18 


clump 


443 


-0.15 


0.079 


19 


-0.21 


clump 


456 


-0.31 


0.206 


27 


-0.21 


RGB 






NGC 


6134 






404 


+0.12 


0.145 


22 


+0.11 


clump 


929 


+0.20 


0.224 


23 


+0.24 


clump 


875 


+0.10 


0.253 


29 


+0.16 


clump 


428 


+0.19 


0.049 


11 


+0.22 


clump 


421 


+0.09 


0.148 


16 


+0.11 


clump 


527 


+0.09 


0.218 


21 


+0.05 


clump 






IC 4651 






27 


+0.16 


0.200 


25 
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clump 


76 


+0.07 


0.141 


19 


+0.10 


clump 


72 


+0.18 


0.172 


21 


+0.13 


RGB 


56 


-0.68 


0.238 


24 


-0.34 


RGBtip 


146 


+0.18 


0.198 


23 


+0.10 


clump 



6.1. Comparison with other determinations 

None of these clusters has a previous metallicity mea- 
sure based on high resolution spectroscopy, but they have 
been the subject of many studies; we present in Table 
|S1 literature metal abundance, based on low resolution 
spectroscopy, or photometric metallicity indicators (in 
DDO, Washington and Stromgren filters), or comparison 



of observed CMDs with theoretical isochrones/tracks. The 
reader is referred to the original papers for detailed expla- 
nations, and we give here only a few notes on some of the 
works. 

Twarog et al. (1997) tried to derive in a homogeneous 
way the properties of a large sample of OC's, and re- 
examined literature data to find distances, reddenings and 
metallicities; for these they selected two methods, DDO 
photometry and the low resolution spectroscopy of Friel 
& Janes (1993), putting the two systems on the same scale. 
Values cited in Table |S1 come from their Tables 1 and 2. 

Friel & Janes (1993) collected low resolution spectra 
of giants of a quite large sample of OC's; for an update, 
see Friel et al. (2002). 

Marconi et al. (1997) used the synthetic CMD tech- 
nique to determine at the same time distance, reddening, 
age, and approximate metallicity of NGC 2506; they em- 
ployed three different sets of evolutionary tracks (Padova, 
Geneva, and FRANEC) finding that, for Z = 0.008 and 
0.01 ([Fe/H] ~ —0.4 and —0.3), they were able to repro- 
duce the observed CMD, although the metal-poor solu- 
tions were preferred. A similar method was employed by 
Meibom et al. (2002) for IC 4651, but they only considered 
the Yale isochrones. 

Finally, Anthony-Twarog & Twarog (2000) give two 
alternative solutions for IC 4651, based on different rela- 
tions for the intrinsic Stromgren colours. 

When we compare our findings with literature values, 
we find i) that our abundance for NGC 2506 is generally 
much higher, ii) that NGC 6134 has strongly discrepant 
determinations, and iii) that IC 4651 is in much better 
agreement with past works. We emphasize, though, that 
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fine abundance analysis of higfi resolution, high signal to 
noise spectra is the most precise method to measure the 
elemental composition. Given also the tests done on our 
temperature and gravity scale, we feel confident about the 
derived metallicities. 

Finally, note that, adopting our metallicity ([Fe/H] = 
+0.15), and the age derived by Bruntt et al. (1999: 0.69 ± 
0.10 Gyr), NGG 6134 appears almost a twin of the Hyades 
for which Ferryman et al. (1998) give [Fe/H]=+0.14 ± 
0.05, age 0.625 ± 0.05 Gyr, and distance modulus (m- 
M)o = 3.33 ± 0.01. The similarity (although the clump 
stars in NGG 6134 are more numerous) is confirmed by 
the absolute magnitude of the clump stars, which span a 
similar range (My ^ 0.2 - 0.5) in the two OG's. 

7. Summary 

We derived precise metallicities for 3 open clusters (NGG 
2506, NGG 6135 and IG 4651) from high resolution spec- 
troscopy, using both the traditional line analysis and ex- 
tensive comparison with synthetic spectra. Our adopted 
temperature scale from excitation equilibrium gives con- 
sistent values of reddenings in very good agreement with 
previous, independent estimates of E{B — V) for the three 
clusters. This finding strongly supports the adopted tem- 
perature scale and, in turn, the derived metallicity scale. 
The nice agreement between spectroscopic and evolution- 
ary gravities also indicates the goodness of the adopted 
temperatures and leaves little space for the effect of pos- 
sible departures from the LTE assumption. Our approach 
is then well suited to derive metal abundances of stellar 
populations with a [Fe/H] ratio around solar. 
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